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Magnetic fields in jets: ordered or disordered? 



R.A. Laing 1 , J.R. Canvin 2 , and A.H. Bridle 3 

1 European Southern Observatory, Karl-Schwarzschild-StraBe 2, 85748 Garching, Germany 

2 School of Physics, University of Sydney, A28, Sydney, NSW 2006, Australia 

3 National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, USA 



Received; accepted; published online 



Abstract. The question of the degree of order in the magnetic fields of relativistic jets is important to any understanding 
of their production. Both vector-ordered (e.g. helical) and disordered, but anisotropic fields can produce the high observed 
degrees of polarization. We outline our models of jets in FRI radio galaxies as decelerating relativistic flows. We then present 
theoretical calculations of the synchrotron emission from different field configurations and compare them with observed 
emission from FRI jets. We show that large-scale helical fields (with significant poloidal and toroidal components) are 
inconsistent with observations. The combination of an ordered toroidal and disordered poloidal component is consistent with 
our data, as is an entirely disordered field. lets must also contain small, but significant amounts of radial field. 
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1. Introduction 

This paper addresses two questions: 

1. Is the magnetic field in extragalactic radio jets vector- 
ordered or does it have many reversals? 

2. What is its three-dimensional structure? 

Observations of polarized synchrotron emission at frequen- 
cies where Faraday rotation is negligible can be used to in- 
fer a two-dimensional projection of the field structure on the 
plane of the sky (we refer to the apparent magnetic field di- 
rection, which is the perpendicular to the observed E-vector 
position angle). This emission is an emissivity-weighted in- 
tegral through the jet, affected by relativistic aberration if the 
flow speed is comparable with c but independent of reversals 
in the field. Measurement of Faraday rotation can potentially 
be used to determine the vector field component along the 
line of sight within a jet, but the thermal plasma responsible 
must be within the jet volume, rather than in front of it, as 
usually appears to be the case (e.g. Laing et al. 2006). 

We have developed models of jets in low -power FR I (Fa- 
naroff & Riley 1974) radio sources as decelerating relativis- 
tic flows. By fitting to deep radio images, we can determine 
the geometry and the distributions of velocity and emissivity. 
Our technique also constrains the three-dimensional struc- 
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ture of the magnetic field. We cannot determine this struc- 
ture uniquely, but we can estimate the ratios of the longitu- 
dinal, toroidal and radial components and eliminate some of 
the simpler proposed field geometries. 

Our fundamental assumption is that jets are intrinsically 
symmetrical and relativistic. Aberration then acts differently 
on radiation from the approaching and receding jets, so their 
observed synchrotron images are two-dimensional projec- 
tions of the field structure viewed from different directions. 
From these, subject to some additional assumptions, we can 
infer the three-dimensional structure of the field. One key as- 
sumption is axisymmetry: if this holds, we can also use the 
symmetry of the transverse brightness and polarization pro- 
files to distinguish (at least in some cases) between vector- 
ordered and disordered fields. 

2. Jet models 

We model FRI jets as intrinsically symmetrical, axisymmet- 
ric, relativistic, stationary flows, in which the magnetic fields 
are assumed to be disordered, but anisotropic (see Section 
for a detailed discussion of this assumption). We adopt sim- 
ple, parameterized functional forms for the geometry and the 
spatial variations of velocity (allowing both deceleration and 
transverse gradients), emissivity and field-component ratios. 
We then optimize the model parameters by fitting to deep 
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Fig. 1. A comparison of VLA data and model for NGC 315 
(Canvin et al. 2005). Vectors with lengths proportional to the 
degree of polarization, p, and directions along the apparent 
magnetic field are superimposed on grey-scales of total in- 
tensity. The resolution is 2.35 arcsec. The polarization and 
angular scales are indicated by the labelled bars in the lower 
panel, (a) model; (b) data. 

VLA images in /, Q and U. The model brightness distribu- 
tions are derived by integration along the line of sight, in- 
cluding the effects of anisotropy in the rest-frame emission, 
aberration and beaming. Note that modelling of the linear po- 
larization is essential to break the degeneracy between angle 
to the line of sight and velocity. Details of the models are 
given by Laing & Bridle (2002), Canvin & Laing (2004) and 
Canvin et al. (2005). A comparison between our best-fitting 
model and VLA data for the first arcminute of the jets in 
NGC 315 is shown in Fig.[T](Canvin et al. 2005). 

3. Field configurations 

The models we have discussed so far fit the observations very 
well, but make a specific assumption about the structure of 
the field: namely that it is disordered on small scales, but 
anisotropic. An alternative class of field structures, proposed 
on theoretical grounds, has a vector-ordered, quasi-helical 




Fig. 2. Sketches of the three field configurations discussed in 
the text. Left: vector-ordered helix. Middle: perpendicular- 
field spine and longitudinal-field shear layer. Right: two- 
dimensional field sheets wrapped around the jet axis. 

structure. In this section, we consider how to differentiate ob- 
servationally between these structures. 

The generic polarization structure which we seek to 
model has an apparent magnetic field which is transverse on 
the axis of the jet but longitudinal at its edges (Fig.Q. Three 
field structures which can all produce this effect are sketched 
in Fig. |3 The first is an ordered, helical field. For simplic- 
ity, we take this to have a constant pitch angle, but qualita- 
tively similar results are obtained from generically similar, 
but more complex configurations. The second has a spine of 
two-dimensional field sheets with equal components in direc- 
tions orthogonal to the jet axis but no longitudinal component 
surrounded by a longitudinal-field shear layer (Laing 1980, 
1993). The third has two-dimensional field sheets wrapped 
around the axis, giving equal toroidal and longitudinal, but 
no radial component (model B of Laing 1980). 

Fig-EJiUustrates the brightness and polarization structure 
produced by a helical field in a jet. In general, helical fields 
produce asymmetric transverse brightness and polarization 
profiles (Laing 1981). The profiles are symmetrical only if 
the field is purely toroidal or the jet is at 90° to the line of 
sight in the rest frame of the emitting material. The condition 
for the approaching jet to appear side-on in the rest frame is 
(3 = cosf? 1 . The counter-jet can never appear edge-on in its 
rest frame unless = and 8 = 90°. Fig |4] shows average 
transverse profiles of total intensity and degree of polarization 
for the jets in NGC 315 (Canvin et al. 2005). Our modelling 

1 (3 = cos 6 is also the condition for maximum Doppler boost, so 
there is a selection effect in favour of observing symmetrical trans- 
verse profiles in blazar jets, even if their fields are helical. 
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Fig. 3. Total intensity and linear polarization from a (non- 
relativistic) jet containing a helical magnetic field of pitch 
angle 45° at angles to the line of sight of 6 = 45° and 90°. 
Left: sketch showing the projection of field lines on the plane 
of the sky. Full line, 9 = 45°; dotted line, 9 = 90°. Mid- 
dle and right: grey-scales of total intensity with superposed 
vectors whose lengths are proportional to the degree of po- 
larization and directions along the apparent magnetic field. 
Middle: 9 = 45°; right: 9 = 90°. 



shows that the angle to the line of sight is « 38° and that the 
field must have a a mixture of longitudinal and toroidal com- 
ponents. The profiles are extremely symmetrical, especially 
in the counter-jet, so we can rule out an ordered helical field 
at least in this source. The toroidal field component could be 
ordered, for example if the longitudinal component has many 
reversals. More generally, a symmetrical transverse profile re- 
quires a symmetrical distribution of field directions in the rest 
frame of the jet flow. 

The other two field geometries sketched in Fig.[2]can be 
distinguished by the relative positions of the transition be- 
tween longitudinal and transverse apparent field on the jet 
axis in the main and counter-jets, provided that the flow is 
relativistic and decelerating. In the case of a perpendicular- 
field spine surrounded by a longitudinal-field shear layer, the 
apparent field transition always occurs further from the nu- 
cleus in the counter-jet and may not be seen in the main jet. 
For field sheets wrapped around the jet axis, giving longitudi- 
nal and toroidal components only, the converse is true. This is 
illustrated in Fig. [5] It has been known for some time that the 
brighter jet in an FR I source tends to show a transition from 
longitudinal apparent field close to the nucleus to transverse 
apparent field further out (Bridle 1984). Our deeper obser- 
vations indicate that the counter-jets either show transverse 
apparent fields on-axis over their entire lengths, or have a 
transition much closer to the nucleus (Laing & Bridle 2002; 
Canvin & Laing 2004; Canvin et al. 2005; see also Hardcas- 
tle et al. 1997). We conclude that a first approximation to the 
field configuration in FR I jets is a mixture of toroidal and 
longitudinal components. 

Detailed modelling (in which the ratios between the field 
components are allowed to vary as functions of position) con- 
firms this general conclusion, but adds detail. Our conclu- 
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Fig. 4. Transverse profiles of total intensity and degree of po- 
larization, p, for the jets of NGC315 (Canvin et al. 2005). 
Full lines: VLA data; dashed lines: model, (a) /, main jet; (b) 
/, counter-jet; (c) p, main jet; (d) p, counter-jet. 

sions for the four sources we have studied in detail are as 
follows: 

- Fields on kpc scales in FR I radio jets are not vector- 
ordered helices. Nor should they be: for a conical jet 
with velocity j3c and bulk Lorentz factor T, the fluxes 
for the longitudinal and transverse field components are 
oc r~ 2 and oc (r/3r) _1 , respectively, where r is the dis- 
tance from the nucleus. The longitudinal flux must be 
very small at large distances from the nucleus, otherwise 
the energy density in the magnetic field close to the ac- 
cretion disk or black hole would be too high (Begelman, 
Blandford & Rees 1984). 

- The field is primarily toroidal and longitudinal, with a 
smaller radial component in some objects. 

- The toroidal component could be ordered, provided that 
the longitudinal component has many reversals. 



References 



References 



(a) Perpendicular-field spine + 
longitudinal-field shear layer 



Field transition 



(b) Toroidal + longitudinal field 




Fig. 6. Panels (a) - (c): grey-scales of the fractional magnetic- 
field components in NGC 315 from Canvin et al. (2005). (a) 
Radial; (b) toroidal; (c) longitudinal, (d) Sketch of the domi- 
nant field components at different locations in the jet, as de- 
duced from the model fits. 
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Fig. 5. Location of the transition between transverse and lon- 
gitudinal apparent magnetic field on the jet axis. The images 
are of model decelerating, relativistic jets (the approaching 
jet is on the right). Vectors representing the degree of po- 
larization and the apparent magnetic field direction are su- 
perposed on a grey-scale of total intensity. The location of 
the field transition is marked, (a) Perpendicular-field spine 
and longitudinal-field shear layer, (b) Field sheets wrapped 
around the jet axis, so that the toroidal and longitudinal com- 
ponents are equal. 

- The longitudinal/toroidal field ratio decreases with dis- 
tance from the nucleus, qualitatively as expected in an 
expanding, decelerating flow. 

- The evolution of the field-component ratios is not, how- 
ever, consistent with flux freezing in a laminar velocity 
field of the type we infer, even if we include the effects of 
velocity shear (Laing & Bridle 2004). 

The field-component structure deduced from our model of 
NGC 315 is shown in Fig.© 

4. Further work 

In principle, our modelling technique should work for jets in 
powerful (FRII) sources on kpc scales and for pc-scale jets. 
There are three main reasons why this is difficult: 

- counter-jets are faint and in some cases the jet/counter-jet 
ratios are very large (presumably the flow is faster); 

- the jets are narrow and 

- more powerful jets show significant sub-structure, over 
which we need to average in order to make a robust com- 
parison of the main and counter-jets. 

The only FR II jets which have been well enough resolved to 
give any idea of their 3D field structures are those in 3C 353 



(Swain, Bridle & Baum 1998). There, a roughly equal mix 
of longitudinal and toroidal field also gives the best fit to the 
polarization data, although the jet velocity is not well con- 
strained. A slightly different approach is feasible in micro- 
quasars, where components from a single ejection event can 
often be identified on both sides of the nucleus. Our technique 
could then be used to compare the brightness and polarization 
of the jet and counter-jet components directly, with the added 
advantage of known proper motions. 

For these new applications, we will need sensitivities of 
0.1-1 /iJy rms and > 10 5 : 1 dynamic range at frequencies 
chosen to be able to resolve and remove Faraday rotation. 
The resolutions required range from 0. 1 arcsec for the nearest 
FRII sources to < 0.1 milliarcsec for parsec-scale jets. We 
look forward to EVLA, eMERLIN and broad-band VLBI. 
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